Heat fluctuations in Brownian transducers 
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Heat fluctuation probability distribution function in Brownian transducers operating between two 
heat reservoirs is studied. We find, both analytically and numerically, that the recently proposed 
Fluctuation Theorem for Heat Exchange [C. Jarzynski and D. K. Wojcik, Phys. Rev. Lett. 92, 
230602 (2004)] has to be modified when the coupling mechanism between both baths is considered. 
We also extend such relation when external work is present. 
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Noncquilibrium systems are receiving much attention 
from a theoretical point of view through the derivation 
of the so called fluctuation theorems. The theoretical ap- 
proach is based on microscopic reversibility and elegant 
analytical properties for the probability distribution of 
entropy generated are derived From these rigorous 
results, corollaries such as a statistical derivation of the 
Second Law can be achieved. Moreover, they hold for sys- 
tems arbitrarily far away from equilibrium and are not 
restricted to the linear regime. There are different fluc- 
tuation relations depending on the dynamics they apply 
to, the magnitudes they relate or the state of the sys- 
tem they refer to. Amongst the most relevant one flnds 
the Gallavoti-Cohen fluctuation theorem the Jarzyn- 
ski equality 0, the formalism for steady-state thermo- 
dynamics an extension of the fluctuation theorem 
[3 and an integral fluctuation theorem 0. Apart from 
their intrinsic value to theoretical physics, a vast number 
of experimental applications based on such results have 
been developed 0| . The benchmark of these theorems are 
nanosystems such as molecular motors Nearly all of 
such relations focus on work fluctuations 1^ ^| . In 
contrast, not much attention has been paid to heat fluc- 
tuations. There are only a few contributions to the topic 
in which only one thermal bath is considered. 
Recently, Jarzynski and Wojcik [ij derived a fluctuation 
theorem for heat exchange (denoted as XFT) between to 
systems initially prepared at different temperatures Ti 
and T2^ then placed in thermal contact with one another 
for a certain time and, flnally, separated again. The the- 
orem states that the probability distribution p{Q) of the 
heat exchange Q satisfies 
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where A/3 = — 1/^27 and fc^ = 1 is taken for sim- 
plicity from now on. Both systems must be prepared 
in equilibrium and then placed in thermal contact with 
one another, for an arbitrary period of time tg, during 
which a net heat Q is transferred. This is a very in- 
teresting result because, in the first place, it focus on 
heat fluctuations and, secondly, because of its universal 
character: the prediction depends only on the two tem- 
peratures and not on any characteristic of the systems. 




Figure 1: Ratchet, pawl and spring Brownian motor. 



Nevertheless, the theorem was derived for the presence 
of heat exchange only, without considering other sources 
of energy such as external work or the energy involved 
in the connecting mechanisms. Then, it is supposed that 
the heat lost by one system is exactly compensated the 
amount energy gained by the other. However, any two 
bodies put in contact need a mechanism that connects 
them and therefore an interaction term should be con- 
sidered and its relevance and effects studied. We may 
also notice that when work comes into play, the transfer 
of heat must be revised because, according the Second 
Law, both baths interchange different amounts of heat. 

The purpose of this work is to study the predictions 
and applicability of the XFT in specific models, such as 
Brownian motors. For a simple mechanical system that 
allows heat exchange between two baths, we show that 
the XFT has to be modified and we present an analytical 
calculation of such modification, which does not depend 
on the details of the connecting mechanism. We also 
propose an extension for the case in which external work 
is present. 

Brownian motors are a set of particularly peculiar ma- 
chines that make use of thermal fluctuations of the en- 
vironment they are immersed in to perform useful work 
[Tsj l. There are many different models being the Feyn- 
man ratchet and pawl device (l^ the paradigm. In gen- 
eral, Feynman-like Brownian motors are built from two 
subsystems immersed respectively in two thermal baths 
at different temperatures and connected by some mech- 
anism. These engines, such as the motor in Ref. JJ], fit 
perfectly the study of heat fiuctuations from the point 
of view of the setting of the XFT. Let's consider the 
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transducer in Fig. ^ It has two degrees of freedom, x 
and y, at different temperatures Ti < T2 and connected 
through a spring. Note that in the colder bath there is 
a saw-toothed wheel which acts as a ratchet potential. 
The Langevin equations of motion of this device in the 
overdamped limit, when setting the friction equal to one. 



i = -d^V,ix,y)~V;{x)+T + Ci{t), (2) 

y = -dyV,ix,y)+^2{t), (3) 

where Vc{x,y) — (fc/2)(x — y)^ is a harmonic coupling, 
Vr{x) is the ratchet potential in p7l | and r is an external 
load. ^i{t) and ^2(i) are independent Gaussian white 
noises of correlation {^i{t)^j{t')) = 2TiSijS{t — t'). 

Heat will flow from one bath to the other through the 
spring. We call Q the heat outcoming from the reservoir 
at T2 ■ Note that due to the internal energy stored in the 
spring, in the ratchet potential and due to the mechanical 
work, Q is not equal to the heat dissipated in the thermal 
bath at Ti. From Sekimoto's energetics scheme |0| one 
can obtain Q as 
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Before a formal approach to the analytic properties of 
this central quantity, let's perform a numerical study of 
it. A numerical computation of the stochastic heat Q 
is easily done by simulating the dynamics of the motor 
and performing the above integral numerically, once the 
steady state has been reached. Then we can get a col- 
lection of values of the total heat Q transferred during a 
fixed time interval ^o, with which histograms can be built. 
This is not exactly the situation described in the deriva- 
tion of the XFT. According to Ref. 0] , both subsystems 
are initially in equilibrium, then connected for a time pe- 
riod to and, finally, separated again. However, in our 
simulations and calculations the transient regime hardly 
contributes compared to the steady state state. Prelimi- 
nary simulations in the nonlinear Brownian transducer at 
nonzero external work conditions show that, in the limit 
kto >> 1, heat histograms are very well fitted by Gaus- 
sian distributions. For smaller kto distributions deviate 
from Gaussianity (see Fig. Several numerical explo- 
rations performed with different values of the parameters 
involved, measuring the kurtosis and the skewness of the 
histograms and checking the tails of the distributions re- 
veal a wide range of the parameter space in which the 
Gaussian approximation is justified. Therefore, in an ap- 
propriate and also physical limit, we can write 
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and, accordingly, the XFT ^ can be expressed as. 
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Figure 2: Histograms of the heat transfer Q for the model 
Q-lO at nonzero external torque for a Gaussian and a non- 
Gaussian case. The parameters are t = 5, Vq = 1, d = 12, 
to = 1, 7i = 1 and T2 = 1.5. The continuous lines are 
Gaussian fits. 



This is a very interesting situation because, if we were 
able to calculate the quantities (Q) and = (AQ^) for 
any specific model, then we could test the XFT prediction 
on a non-ideal case. This cannot be done analytically for 
the nonlinear model but it is possible if we sim- 

plify it by taking Vr — 0. In particular, for this passive 
Brownian transducer and in the absence of external load 
(r — 0) we find 



(Q)=^o^(T2-Ti), 
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which is indeed Fourier's law for the thermal conductivity 
(see Refs. Ij9 |.l20l ^ and, in the limit kto >> Ij the second 
moment is |2l| 



(8) 



Note that the mean value and variance are extensive in 
to and they are also linear in k. The above results yield 
to 



F(0) = ln^i±^ = A/3Q(l-7), 



where 
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This is one of the main results of this work. The ex- 
pression does not depend on any detail of the coupling 
mechanism. We obtain that the XFT holds for small 
A/3 {T2 —Ti) but an important correction of order Af3'^ 
appears. In Fig. Owe plot the XFT the new pre- 
diction for the passive Brownian transducer ||5J) and re- 
sults from numerical simulations, exploring the T2 depen- 
dence. The symbols correspond to numerical data and 
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Figure 3: Prediction of the XFT (equation and the cor- 
rection found due to the coupling mechanism (equation ||5J) 
in the Gaussian approximation. The symbols are obtained 
from numerical data through different procedures (see text). 
The device considered is the passive Brownian transducer at 
r = and Ti = 1. 



Figure 4: Numerical simulations (dots) and analytical predic- 
tion (line) of equation for the relation of heat fluctuations 
versus the mean value of the work in the Brownian motor 
model defined in equations The agreement is very 

good, even for this nonlinear system. The data are obtained 
by direct measurement of the slope of Y{Q). The parameters 
used are Ti ^1,T2 = 1.5, d = 12, Vo = 1, to ^ 1 and k = 100. 



deserve a careful explanation. The circles are obtained at 
kto — 100. For white circles, Y{Q) has been obtained by 
direct analysis of the probabihty distribution functions of 
heat. Around T2 > 2, negative values {Q < 0) are rarer 
events and therefore the tails are very difficult to observe 
directly. Black circles are obtained assuming Gaussian 
behavior of the tails, though it is important to note that 
the distributions are still checked to be fitted reasonably 
well by a Gaussian. The higher the temperature T2, the 
worse the Gaussian supposition is and so, for T2 ^ 6, 
we cannot conclude anything about Y{Q). Finally, dia- 
monds correspond to data at kto — 10 from which Y{Q) 
is directly measured from the histograms. This case is 
also very instructive because the distributions observed 
are non-Gaussian (see Fig. The prediction ^ has 
also been checked satisfactorily (data not shown) in the 
nonlinear Brownian motor for kto = 100, Vg = 1 and 
d = 12. Then our simulations indicate that, despite the 
rapidly increasing difficulty to test the theorem, impor- 
tant and systematic deviations do appear. 

The 7 factor modifying the XFT is a signature of the 
important role of the coupling between both systems. 
This points out the applicability of the theorem for large 
A/3, which strongly depends on the approximation made 
in the XFT derivation. It consists on neglecting the inter- 
action term coupling the two bodies by assuming that 
the energy involved in the coupling mechanism is much 
smaller than the typical energy change in both systems. 
When comparing in our model the typical energy of the 
coupling mechanism (the mean value of the potential en- 
ergy of the spring (V^) = {T2 + Ti)/4), and the typical en- 
ergy change in every subsystem (the mean heat released 
{Q}), for the parameter values k = 100, to — 1, T2 — 5 
and Ti = 1, we find (Vc)/(Q) - 0.0075. Therefore, we 
must remark here that, albeit the interaction energy can 
be negligible, its consequences are not. 



The second point we study is the effect of an external 
work into the system. We want to explore the possibility 
of an extension of the relation (jSj) for this case. Taking 
T 7^ in (PI, discarding the nonlinear ratchet potential 
and proceeding as before for these type of calculations, 
the quantity (Q) is 
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where {Q)c is the mean heat conducted and (W) = 
is the mean work. Remember that (Q) is the mean 
heat released by the heat bath a T2, so we are studying 
the fluctuations of this quantity, while heat exchanged at 
bath Ti is different. The Fourier heat is conducted to 
the cold bath but also the hotter bath receives heat from 
the colder due to the external work. In fact, each bath 
dissipates half of the total work. It is worth remarking 
that the sign of (Q) can be reversed for (W) > 2{Q)c, 
reversing the heat flux, now from the cold bath to the 
hot one. The calculus of the variance is more involved 
but we find 
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Using these results and defining the ratio of the two 
(model dependent) energies involved, TZ = {W)/{Q)c, 
the generalized relation for heat fluctuations is 
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The above extension depends on the mechanisms in- 
volved. Nevertheless, although, expression l(T^ has been 
derived for a linear model, we can apply it to the nonlin- 
ear case obtaining a very good agreement, as it is shown 
in Fig. 0] This means that we have found the terms that 
gather the most relevant features and which work even 
for general nonlinear devices. Notice that the torque used 
is in general beyond the stall torque of the motor perfor- 
mance. This is so because for very small loads, thus the 
ones that this motor is able to lift, it is very difficult to 
observe changes in the distributions of Q. We must also 
stress that our analytical prediction does not involve any 
adjustable parameter and, as a consequence, it could be 
confronted against other type of conducting and working 
devices. 

One can derive similar results if the heat interchanged 
by the cold bath at Ti is considered instead of the heat 
transferred from the bath at T2. In this case, we have 
found (not shown here) that the relation is unchanged 
when we deal with the fluctuations of the energy dissi- 
pated in the cold bath. Nevertheless for the t 7^ case, 
the expressions vary but one can find the corresponding 
relations following the same type of calculations. With 
respect to the conditions of our approach, we stress that 
it could be possible to obtain analytical expressions for 
the statistical moments of Q in the transient regime, 
after putting both systems initially in contact. Tran- 
sient corrections of order e~^^° appear which, in the limit 
kto » 1, can be discarded. Thus the transient regime 
is negligible compared to the steady state contribution. 
Therefore, all the calculations in this letter are done in 
the steady state and in the long limit (or big coupling) . 
This is, indeed, of great advantage because it makes it 
possible to derive analytically the most dominant contri- 
butions of the first and second moments. What is more. 



in such limit, although p{Q) is not rigorously Gaussian 
because Q{t) as defined in is a nonlinear functional 
of a Gaussian Orstein-Uhlenbcck process, our main re- 
sults are dominated by the Gaussian-like property of the 
distribution. As a byproduct we have shown in a linear 
model that heat fluctuations relative to the mean value 
cr^/((5) are system independent. It would be worthy to 
explore this result for other nonlinear models. 

This work is the first application and test of the XFT 
U to a non-ideal system, in the sense that the effects of 
the system-environment coupling energy cannot be ne- 
glected. We have checked the sensitivity of the hypothesis 
of small interaction term in Ref. T| for simple Brownian 
devices. Quite surprisingly the XFT works better when 
the heat conducted is of the order of the energy stored in 
the coupling device but not when it is larger. It is worth 
emphasizing that any attempt to write a model that con- 
sists on two bodies interchanging heat needs a connection 
and, therefore, the observations mentioned above are en- 
countered. Hence, one can neglect the energy stored in 
the connecting mechanism but its effect has relevant con- 
sequences. However, in order to understand in more de- 
tail the role of heat, work and coupling energy, it would 
be very interesting to address these questions from a more 
general theoretical point of view. The applicability of 
such results in theoretical models and in experiments is 
of great importance for discovering and understanding 
noncquilibrium statistical mechanics principles. 
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